Astronomy & Astrophysics manuscript no. 48 1 1 main
(DOI: will be inserted by hand later)
March 3, 2009
GRB 030406 - an extremely hard burst
outside of the INTEGRAL* field of view
R. Marcinkowski 1 - 2 , M. Denis 1 , T. Bulik 3 ' 4 , P. Goldoni 5 - 6 , Ph. Laurent 5 ' 6 , and A. Rau 7
O
O
(N
u
Oh
<
O
>
m
o
O
Or
i '
O'
a:
\
Space Research Center, Bartycka 18a, Warsaw, Poland
IPJ, 05-400 Swierk/Otwock, Poland
Astronomical Observatory, Warsaw University, Aleje Ujazdowskie 4, 00478 Warsaw, Poland
Nicolaus Copernicus Astronomical Center, Bartycka 18, 00716 Warsaw, Poland
SAp CEA, 91911 Gif-sur-Yvette, France
UMR 7164, 1 1 Place M. Berthelot, 75231 Paris, France
Max-Planck Institute for extraterrestrial Physics, Giessenbachstrasse, 85748 Garching, Germany
Received date / Accepted date
Abstract. Using the IBIS Compton mode, the INTEGRAL satellite is able to detect and localize bright and hard GRBs,
which happen outside of the nominal INTEGRAL field of view. We have developed a method of analyzing such INTEGRAL
data to obtain the burst location and spectra. We present the results for the case of GRB 030406. The burst is localized with the
Compton events, and the location is consistent with the previous Interplanetary Network position. A spectral analysis is possible
by detailed modeling of the detector response for such a far off-axis source with an offset of 36.9 °. The average spectrum of
the burst is extremely hard: the photon index above 400 keV is -1.7, with no evidence of a break up to 1.1 MeV at a 90%
confidence level.
Key words, gamma rays: burst
1. Introduction
The INTEGRAL satellite detects gamma-ray bursts (GRBs)
in two different ways: for a small number of events
that f all in the field o f view (FoV) of the imager
IBIS (Uberti nTeTail 12003) and of the spectrometer SPI
Vedrenne et al. 2003), INTEGRAL provides accurate posi-
tions (~2 arcmi n) for rapid ground- an d space-based follow-
up observations (Mereg hetti et al.l2003l) . A significantly larger
number of GRBs occurs outside of the FoV of the two
instruments. These bursts can be monitored by the SPI
anti-coincidenc e system: SPI- ACS dvon Kienlin et alJ |2003;
iRau et al.ll2005l) . Up to now, the localization of the bursts that
happened outside of t he field of vie w was only possible with
the aid of the 3 rd IPN faurlevl 19971) .
In this paper we show that for some of these bursts it is
also possible to perform a more detailed localization analysis
using the Compton mode of IBIS, provided that the burst is
sufficiently strong and spectrally hard.
Send offprint requests to: R. Marcinkowski, e-mail:
radek@cbk . waw . pi
* INTEGRAL is an ESA project with instruments and science
data center funded by the ESA members (especially the PI coun-
tries: Denmark, France, Germany, Italy, Switzerland, Spain), Czech
Republic, and Poland, and with the participation of Russia and USA
We present the capabilities of the Compton mode using
an example of GRB 030406, a burst that was detected out-
side of the INTEGRAL field of view. GRB 030406 has been
detected by SPI-ACS on-board INTEG RAL as wel l as by
Ulysses, Konus, and Mars Odyssey (Hurlev et al. 2003). It was
a long burst lasting ~ 65 seconds. Ulysses reported the fluence
of 1.3 X 10~ 5 ergcrrT 2 (25-100 keV). The observations by the
three satellites yielded an estimate of the position centered on
RA(2000) = 19'' l m 43.0 s , DEC(2000) = -68° 4' 39.4", with
the 3cr error-box of 77 square arcminutes.
In section |2 we describe the Compton mode of IBIS
INTEGRAL and demonstrate its localization capabilities. In
section 12.21 we present the spectral analysis of GRB 030406
while section|3]contains the discussion.
2. Observations and data analysis
The IBIS telescope is an imaging instrum ent on-board the
INTEGRAL satellite with a coded mask ( Goldwurm etalJ
20031). There are two dete ction layers i n the detector plan e
ISGRI JLebrun et alJl2003l) and PICsIT JLabanti et al.ll2003h .
ISGRI is an array (128 x 128) of pixels made of semiconduc-
tive CdTe, sensitive to photons between 15keV and ~1 MeV.
ISGRI works in photon-by-photon mode. PICsIT having the
same detection area as ISGRI, is an array of 64 x 64 Csl scintil-
lators, sensitive between ~ 1 70 ke V and 1 5 MeV located 94 mm
R.Marcinkowski, M.Denis, T.Bulik: GRB 030406
below ISGRI. PICsIT (usually) makes light curves of the whole
detector by making histograms of all events (so called Spectral
Timing Mode). ISGRI and PICsIT can act as a Compton tele-
scope, registering photons that are scattered in one and ab-
sorbed in the other detector. The coincidence time window is
a parameter programmable on-board and it was set to ~ 4 mi-
croseconds at the time of GRB 030406. The Compton mode is
sensitive between 200 ke V and ~5 Me V.
There is a 3 meter long collimator on top of the detector
unit. The walls of the collimator are made of lead and act as a
shield to photons with energies up to » 200 ke V. For the geo-
metrical reasons the optical depth of the shield is smaller for
photons arriving at large angles. Thus hard photons from off-
axis sources can pass through the shield and reach the detector.
In particular off axis GRBs may be detectable in IBIS.
Thanks to the Compton mode data we were able to make
successful, independent burst localization. Including informa-
tion from the ISGRI data we have performed burst spectral
analysis.
2.1. Compton imaging and localization
We have analyzed the Compton mode count rate during the
GRB 030406. The burst profile in the Compton mode fol-
lows clearly the profiles seen by SPI-ACS, ISGRI and PICsIT
(Fig. Q. The gamma ray burst was detected in the Compton
mode at the level of 30cr reaching ~400 ctss -1 in the sec-
ond peak (average background amount to ~100 ctss -1 ). The
Compton mode telemetry did not suffer the data gaps in con-
trast to ISGRI.
The Compton mode provides us with the following infor-
mation:
- energy deposits in ISGRI (Ei) and PICsIT (£ P ),
- position of the detection in ISGRI (xy, yi) and PICsIT (xp, yp),
- timing of the event.
The two positions provide the information about the di-
rection of the scattered photon. The two energies provide the
Compton scatter angle 6c, which cosine, in case of the forward
scatter, is given by:
cos 8c — 1 -
m e <r
E\ E\ + Ep
where m e c 2 is the electron rest energy. Given the angle 0c and
the direction after scattering we obtain a ring in the sky which
contains the possible directions of the primary photon. In case
of a point source all the rings cross at the location of the source.
For a given observation a collection of such rings provides a
Compton map of the sky S . In real observations the Compton
map S consists of the source signal and the background B.
In the case of GRB 030406 we selected the Compton events
from the begin of the peak up to the end of the tail, which
consist of a 9.75 s time interval, see Fig.^ For these events a
Compton map S y of the sky was calculated using the method
described above with the angular resolution of one degree. In
order to determine the background we used the pre-burst data
divided into 200 intervals lasting 9.75 s each. For each interval
we calculated the all-sky Compton map and we determined the
mean background B,j and its variance Vtj in each pixel.
L K J
"i ' ' r
u SPI
1 ACS
3
o
u
(H
IBIS
ISGRI
Kr
IBIS
Compton mode
140000
120000
100000
80000
1000
_500
X
300
200
100
4000
2000
10 20 30 40 50 60
Time after trigger [s]
Fig. 1. Lightcurve of GRB 030406 seen by the INTEGRAL
sub-systems: SPI-ACS, IBIS/ISGRI, IBIS Compton and
IBIS/PICsIT The gray areas: precursor, peak and tail corre-
spond to the time regions chosen for spectral analysis. Time =
corresponds to the SPI-ACS trigger time at UTC 2003-04-06,
22:42:03.23.
The map <x /; was obtained by dividing the difference be-
tween the Compton map S u and the background B (/ by the map
of the square root of variance of the noise in each pixel JVij.
We present the map of significance <x (i in Fig. [2] The largest de-
viations from the noise are ~20, and the position of the pixels
with the strongest deviation is consistent with the IPN location.
In order to obtain an estimate of the accuracy of such posi-
tioning we used two approaches: the maximum likelihood and
the Monte-Carlo. In the first approach we assumed that the
probability density that the true position lies in a given pixel
i, j is proportional to the likelihood function:
IBIS
PICSIT SPTI
£ u oc exp (cr, 7 ) .
After normalizing such probability density we found the region
containing 68% of the probability corresponding to the ~ 3.5°
radius.
In the se cond approach we used the IBIS mass model
(Laurent et a l. 2003) to extract the signal from simulated GRB
R.Marcinkowski, M.Denis, T.Bulik: GRB 030406
5
O
80
60
5^ 40
20
a
|-20
m
Q-40
-60
-80
[°]
10-
o
10
I
-20 -10 10 [°]
60
120
180
RA[°]
240
300
Fig. 2. All-sky Compton map taken during the second peak of the GRB 030406. The inset shows the dashed area around the
GRB 030406 centered on the IPN position.
1
i
i i
yj/o
-
0.8
68% ,
|0.6
1-7
J
0.4
S
-
0.2
— 1 i
i
i i
2 4 6 8 10
Position offset [°]
Fig. 3. The normalized cumulative distribution of the offset of
best fit localizations of simulated bursts obtained in the Monte
Carlo simulation: 68% of the events are localized within ~ 4°
of the actual position.
events inserted into stretches of data with experimental back-
ground. We then performed the identical data analysis proce-
dure and found the localization of such simulated bursts. We
have analyzed 127 simulated bursts, and found the cumula-
tive distribution of the distance between the actual burst po-
sition and the estimate obtained in our analysis. This distribu-
tion is presented in Fig. [3] The Icr region (containing 68% of
the cases), corresponds to the accuracy of ~4°, and the 2cr is
constrained to 6.5°. The two methods lead to consistent results.
2.2. Spectral analysis
In addition to the Compton mode GRB 030406 has been also
detected by ISGRI and PICsIT There are only two energy
channels in the PICsIT data stream therefore we have decided
to model together the Compton and ISGRI data, neglecting the
PICsIT.
We have divided the data into three time intervals marked in
gray in Figure[Oas: precursor, peak and tail. The division of the
main peak of the burst into the peak and tail was a consequence
of the ISGRI data loss during the burst due to telemetry gap.
We first ran Monte Carlo simulations using the mass model
of the INTEGRAL spacecraft to generate the detector response
matrices (DRM) for a source located at the position of the burst.
The DRM matrices were obtained for the ISGRI and Compton
sub-systems. For all of them we have taken into account the
influence of the IBIS VETO system, the energy thresholds in
each ISGRI pixel and PICsIT crystal, the ISGRI noisy pix-
els, the failed ISGRI pixels and PICsIT crystals as well as the
ISGRI and PICsIT energy resolution.
The Compton data were binned into six energy channels
covering the range from 200 to 2500 ke V. The ISGRI data spans
from 24 to 450 keV which partially overlaps the Compton en-
ergy range and were binned into seven channels. The analysis
of the pre-burst data with the correction for the varying dead
time effect and exposures provided an estimate of the back-
ground. This background was subtracted from the total count
spectrum in each channel.
We used the standard XSPEC 1 1 .2 JArnaudll996h to fit the
data. For each time region we tried to fit several models: sin-
gle power law, broken power law, and a non physical model
of single black-body. The single power law model never fitted
the combined ISGRI and Compton data. The broken power law
fits the data very well in all analyzed regions: precursor, peak
precursor
R.Marcinkowski, M.Denis, T.Bulik: GRB 030406
peak
teii
%
- ^J
f^-f*U, 1
"^
•*N
i
o o
1
e °
u
o
o
1
T
3
*°
- — —
— -
—
-4 -
_— J
-
Y
w
:
100 1000
channel energy (keV)
100 1000
channel energy (keV)
100 1000
channel energy (keV)
Fig. 4. GRB 030406 spectrum fits to the combined ISGRI and Compton mode data. Upper panels show count spectra, while lower
corresponding vF v spectra. Broken power law model was fitted in each case: to the precursor data (left), to the peak part data
(center) and to the tail (right). See text for details
and tail. In additional fit we have allowed for different normal-
ization of ISGRI and Compton components for peak region,
taking to account two facts related with dead time influence:
- during precursor and tail the count rate grows with factor
about 1 .5 for ISGRI and Compton mode, while during peak,
these factors reach 2.5 for ISGRI and 4 for Compton mode,
- the width of the burst peak, «ls, is much shorter than the
INTEGRAL dead time effect calculation bin, equal 8s.
Except meaningful decreasing of the reduced x 2 we have not
find significant evidence of changing fitted parameters. Results
of all fits are presented in Tabled The broken power law fits
and the residuals are presented in the top panels of Figure|4]
We note that the precursor and the tail region can also be
fitted with black body spectra with temperatures w lOOkeV.
This is due to the fact that at high energies the count rates fall
steeply. The lower panels of Figure 0] show the deconvolved
vF v spectra in each time interval. The precursor and tail spectra
have peak power around 300-500 keV, while the vF v spectrum
in the peak seems to rise all the way into high energies.
The high energy index of yS = -1.7 + 0.3 for peak may
imply the absence of a peak in the source power spectrum vF v
in the range covered by our data, i.e. up to 2.5 MeV. In order
to estimate the lower limit on such a break compatible with the
data we conducted the following analysis. We fit the Compton
data with a broken power law with fixed high energy index:
f(E) = N
*-* "^ *-* break
break
(i)
t^ •> tL break
In this initial fit the break energy was constant
£break=3 000keV. This value was chosen high enough not
to influence the power law part of the fit. For this fit we
obtained x 2 - 0.55 with 4 degrees of freedom. We used the
standard model parameter estimation dLampton et alJ 1 1976)
to calculate the values of x 1 corresponding to 90% and 99%
confidence regions, which are respectively Ax 2 — 8.3 and
Ax 2 = 13.8. Then keeping all the parameters but Ebreak
constant we looked for the lowest value of £break for which
X 2 was equal to the above values. We present the results in
Figure |5] We obtained £break >1 HOkeV at 90% confidence
level and £break >880keV at 99% confidence level.
3. Discussion
We have shown that the IBIS Compton mode is able to de-
tect and localize GRBs outside of the field of view of the
R.Marcinkowski, M.Denis, T.Bulik: GRB 030406
Table 1. Results of the spectrum fits for three time ranges of the GRB 030406. All the errors are at the lcr level.
Part of duration
the GRB [s]
precursor
peak
peak f
tail
7.3
2.81
4.3
00+51 9X)
-0.3
broken power law fit blackbody fit
ft JWjkeV] Ar 2 /dof kT [keV] ,y 2 /dof
-1.5 +
;+0.7 i 7+O.4
1.0 •'-0.3
-2.5 + '^ 1.7 +0 ' 2
-0.5
-0.2
-0.8!?
+(1.7
2
2-8-oJ
490+ 40
390!*
300f*>
27 +70
0.96
1.19
0.49
0.56
-Wjf
140! 2 jj
0.82
0.56
NOTE: ' data fitted with variable normalization for both sub-instruments. See text.
\
99%
-
^
90%
1000
Ebreak [keV]
Fig. 5. Confidence limits on ^cutoff : tne solid line represents the
X 2 as a function of iicutoff using the model of equation Q. We
present the x 2 values corresponding to the 90% and 99% con-
fidence levels.
INTEGRAL telescopes. We have found that the position of
GRB 030406 obtained using the Compton mode technique is
consistent with the IPN localization. The position accuracy has
been estimated using the maximum likelihood and the Monte
Carlo method. The results of both methods are consistent and
the lcr errors are ~ 4°.
Given the position of the burst we computed the DRM for
a source outside of the field of view of the satellite at this par-
ticular location. Using the DRM we analyzed the the Compton
data, along with the ISGRI detection and fitted the spectrum
with a broken power law in the range from 50 ke V to 3 Me V.
The peak spectrum is very hard; in the vF v the low energy
index below 400 keV rises with the index a +3.5 and above
this energy it is still positive » +0.3. We note that this spec-
trum is an average over the time interval of 2.8s around the
maximum, thus the spectrum of the maximum might have even
been harder since the burst spectra usually evolve from hard in
the peaks to soft in their tails. We do not see any evidence that
the vF v spectrum peaks below 1 . 1 MeV. Thus, there is a hint
of existence of bursts with th e peak of vF v spectrum above the
distribution shown bv lMallozzi et alJ Jl995) which truncates at
~ 1 MeV for the bright BATSE bursts.
Such hard spe ctra have already been seen by BATSE -
Pre ece et al.l (12000) present the distribution of the low and high
energy spectral indexes for a sample of 156 bright bursts. Both,
the low and the high energy spectral indexes of GRB 030406
lie in the upper tail s of the distrib utions of the spectral indexes
presented bv lPreece et alJ 1120001) . It is therefore clear that the
spectrum of this particular burst is in clear contradiction with
the synchrotron model of GRBs (Katz 1994; Tavani 1995),
which predicts that there is a strict upper limit on the low en-
7 J T ^J.
ergy spectral index of -j (Preece et al. 1998). The low energy
spectral slope is consistent within the error bars with the jitter
synchrotron model of GRBs (Medvedev 2000). We note that a
detailed study of time resolved spectra (Ghirland a et alJl2003l)
showed that the low energy spectral slopes are large ~ 0.5 - 1
at the rising parts of several bursts.
The value of E p of GRB 030406 - the peak of the vF v spec-
trum lies above 1 . 1 MeV. It should be noted that thi s is on the
high end of the Amati and Ghirlanda relations l l Amati et alJ
20021 iGhirlanda et aljEooH) . If this relation holds and we as-
sume that the peak energy is larger than the 90% lower limit of
1 lOOkeV, then the isotropic radiated energy for GRB 030406 is
larger than (6 - 20) ■ 10 53 erg, where the uncertainty stems from
the inaccuracy of the Amati relation. Hard bursts are therefore
very energetic and may come from a population residing at high
redshifts.
In summary, we confirm that the INTEGRAL in the
Compton mode can detect hard GRBs. The localization accu-
racy is a few degrees. The GRB spectra can be studied in the
range from ~200keV to ~3 MeV depending on the location of
the burst in the instrument coordinates. Further study of such
hard bursts detectable in the INTEGRAL Compton mode may
yield potential candidates of high redshift burst and also test the
limits of validity if the Amati relation. Our rough estimations
of the Compton mode capabilities give 2-3 detections of such
bursts per year. They are consistent with BATSE results and
seem to be confirmed by the other detections (Marcinkowski et
al. in preparation). Finally, we note that analysis of Compton
data in principle allows to determine the polarization. A study
of of the polarization sensitivity of the IBIS Compton mode is
currently under way (Laurent et al. in preparation).
Acknowledgements. This research was supported by the KBN grants
2 P03D 001 25, 1 P03D 003 27 and PBZ-KBN-054/P03/2001.
References
Amati, L., Frontera, F, Tavani, M., et al. 2002, A&A, 390, 81
Arnaud, K. A. 1996, in ASP Conf. Ser. 101 : Astronomical Data
Analysis Software and Systems V, 17-+
Ghirlanda, G, Celotti, A., & Ghisellini, G. 2003, A&A, 406,
879
R.Marcinkowski, M.Denis, T.Bulik: GRB 030406
Ghirlanda, G., Ghisellini, G., & Lazzati, D. 2004, ApJ, 616,
331
Goldwurm, A., David, P., Foschini, L., et al. 2003, A&A, 411,
L223
Hurley, K. 1997, in ESA SP-382: The Transparent Universe,
491-+
Hurley, K., Cline, T, Mitrofanov, I., et al. 2003, "GCN", 2127,
1
Katz, J. I. 1994, ApJ, 432, LI 07
Labanti, C, Di Cocco, G, Ferro, G, et al. 2003, A&A, 411,
L149
Lampton, M., Margon, B., & Bowyer, S. 1976, ApJ, 208, 177
Laurent, P., Limousin, O., Cadolle-Bel, M., et al. 2003, A&A,
411,L185
Lebrun, F, Leray, J. P., Lavocat, P., et al. 2003, A&A, 411,
L141
Mallozzi, R. S., Paciesas, W. S., Pendleton, G N., et al. 1995,
ApJ, 454, 597
Medvedev, M. V. 2000, ApJ, 540, 704
Mereghetti, S., Gotz, D., Borkowski, J., Walter, R., & Pedersen,
H. 2003, A&A, 411, L291
Preece, R. D., Briggs, M. S., Mallozzi, R. S., et al. 1998, ApJ,
506, L23
— . 2000, ApJS, 126, 19
Rau, A., von Kienlin, A., Hurley, K., & Lichti, G G 2005,
A&A, 438, 1175
Tavani, M. 1995, Ap&SS, 231, 181
Ubertini, P., Lebrun, F, Di Cocco, G, et al. 2003, A&A, 411,
L131
Vedrenne, G., Roques, J.-P, Schonfelder, V, et al. 2003, A&A,
411, L63
von Kienlin, A., Beckmann, V., Rau, A., et al. 2003, A&A, 411,
L299